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PLANAR LINKAGES FOLLOWING A PRESCRIBED MOTION 


MATTEO GALLET* *’°, CHRISTOPH KOUTSCHAN*, ZIJIA LP, 
GEORG REGENSBURGERt, JOSEF SCHICHO, AND NELLY VILLAMIZAR 


Abstract. Designing mechanical devices, called linkages, that draw a given 
plane curve has been a topic that interested engineers and mathematicians 
for hundreds of years, and recently also computer scientists. Already in 1876, 
Kempe proposed a procedure for solving the problem in full generality, but his 
constructions tend to be extremely complicated. We provide a novel algorithm 
that produces much simpler linkages, but works only for parametric curves. 
Our approach is to transform the problem into a factorization task over some 
noncommutative algebra. We show how to compute such a factorization, and 
how to use it to construct a linkage tracing a given curve. 


1. Introduction 

Kempe’s Universality Theorem El. stating that any plane algebraic curve can 
be drawn by a mechanical linkage with only rotational joints, surprised his con¬ 
temporaries: his work represented a major breakthrough in a topic investigated by 
mathematicians of the 19th century. In the middle of the last century, it appeared 
in textbooks, for example, the ones by Lebesgue j^Hl and Blaschke and Muller [5]. 
Recently, there was a revived interest in the problem among mathematicians and 
computer scientists as we outline in the following. We would also like to mention [52] 
for a historical overview and further references. 

In modern terms, the procedure proposed by Kempe is a parsing algorithm. It 
takes the defining polynomial of a plane curve as input and realizes arithmetic op¬ 
erations via certain elementary linkages. In this work, we approach the question 
from a different perspective. Instead of a polynomial, we start with a rational 
parametrization of a curve. From this, we obtain a parametrized family of elements 
of the group SE 2 of direct isometries of the plane, namely a motion, whose action 
on a point traces the given curve. In order to realize such a motion by a linkage, we 
decompose it into a series of revolutions. By encoding motions via polynomials over 
a noncommutative algebra, we reduce this task to a factorization problem. Even¬ 
tually, we design a linkage whose rotational joints move according to the previously 
obtained revolutions, yielding a device drawing the desired curve. 

Requiring a rational parametrization, it is clear that our approach is less general 
than Kempe’s, and its subsequent generalizations [Hiiiiniiin], because it cannot 
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be applied to planar curves of positive genus. However, the parametric setup has 
advantages in design processes: in the language of robotics/kinematics, we are not 
just prescribing the position of the end effector but also its orientation, in a similar 
way as in |25| ; moreover, when one has to interpolate prescribed poses of the end 
effector, we can specify a tracing order and have control over its speed. 

Linkages obtained via Kempe’s procedure do not, in general, trace only the curve 
they are designed for. More precisely, the curve is only drawn by a component 
of the configuration space (also called workspace), that is, the set of all positions 
reachable by the linkage. This is due to the fact that already the elementary linkages 
may admit degenerate configurations, allowing the device to flip into an unwanted 
component. Several solutions to this problem have been suggested [Hlelinilisilll]. 
The linkages produced by our method present the same issue, however, this can be 
treated using the techniques presented in mii; see Remark |7.7[ 

Generation of linkages drawing an arbitrary curve may involve many links and 
joints. This was already observed by Kempe mi: “.. this method would not be 
practically useful on account of the complexity of the link work employed .. The 
problem of finding upper bounds for the number of joints needed to trace a given 
planar curve of degree d was first addressed in 13 [S] ; see [5J Section 3.2.3]. Their 
upper bound was later improved by Abbott and Barton to O(d^), who in [T] 

showed that this bound is optimal. In this paper, we provide an algorithm that ap¬ 
plies to rational planar curves: given a rational parametrization with denominator 
of degree d and without real roots, our a lgorithm produces a linkage with 3d + 2 

we are grateful to Hans-Peter Schrocker 


d + I joints (Proposition 


7.10 


links and | 

for discussions leading to this result). If we apply our technique to the parametriza¬ 
tion of an ellipse, we obtain a linkage with 8 links and 10 joints; in Section]^ we 
illustrate our approach and the main ideas of our paper via this example. In con¬ 
trast, an unoptimized linkage returned by Kempe’s procedure has 158 links and 235 
joints; we are grateful to Alexander Kobel, who gave a full implementation |16| of 
Kempe’s procedure, for assisting us in the computation of these numbers. 

In Section we develop a mathematical model for linkages with only rotational 
joints. In particular, we define their configuration space in terms of isometries send¬ 
ing a fixed initial configuration to a reachable one. This differs from the commonly 
used models miniiisiiigEi] which, instead, define configurations in terms of the 
positions of links and joints. 

In Section]^ we recall that one can embed SE 2 as an open subset of a real projec¬ 
tive space; see m- This allows us to introduce a noncommutative algebra K whose 
multiplication corresponds to the group operation in SE 2 , hence mimicking the role 
played by dual quaternions with respect to SE 3 . A polynomial with coefficients in K 
therefore describes a family of direct isometries, which we call a rational motion. 
Consequently, we refer to such polynomials as motion polynomials; they are the 
two-dimensional analog of the motion polynomials introduced in m, and they are 
the key for turning our geometric problem into an algebraic one. In fact, we show 
that a linear motion polynomial represents a motion constrained by a revolute or 
prismatic joint. Hence a factorization of a motion polynomial into linear factors 
yields a decomposition of the corresponding rational motion into simple ones. 

In Section we give necessary and sufficient criteria for the existence of a fac¬ 
torization of a motion polynomial. Not every motion polynomial admits such a 
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factorization. However, the correspondence between motion polynomials and ratio¬ 
nal motions is not one-to-one: for every rational motion, there is a whole equivalence 
class of motion polynomials. Moreover, the motions we are dealing with are spe¬ 
cial, since they admit bounded orbits. Once we restrict to such bounded rational 
motions, we can prove that the equivalence class of every bounded rational motion 
contains a motion polynomial admitting a factorization into linear polynomials of 
revolute type. We provide an algorithm for computing the smallest factorizable 
polynomial in an equivalence class and one of its factorizations; see Theorem |5.15| 
In this context, we would like to refer to the preprint m, where this result is used 
to prove factorizability of motion polynomials giving rational motions in SE 3 . 

At the beginning of Section we construct from a factorization a linkage, in 
the form of an open chain, such that the given rational motion can be realized 
as the relative motion of the last link of the chain with respect to the first one. 
Open chains have high mobility, so we have to constrain our linkage such that it 
performs only the motion we are interested in. The technique we employ, called 
flip procedure, is introduced in Section]^ In the remaining part of Section ITlwe 
exploit the properties of flips and propose our main algorithm; see Theorem |7.6[ 

In Section]^ we address the problem of self-collisions. For arbitrary linkages, this 
is a challenging problem 0 Section 9.3]. Here, for the first time in the paper, we take 
into account how linkages are physically realized. We show that self-collisions can 
be efficiently detected for linkages obtained by our algorithm when links are realized 
by bars. If, instead, we allow links of different shapes, we describe a construction 
showing that it is possible to realize these linkages without collisions. This addresses 
the open problems raised in Open Problem 3.2] and |23[ Section 2.3]. 

A popular formulation of Kempe’s Theorem states that “There is a linkage that 
signs your name”. Inspired by this, in Section]^ we give an example of a linkage 
drawing a calligraphic letter. 

All algorithms described in the paper have been implemented by the second- 
named author in the computer algebra system Mathematica. The source code, an 
expository notebook, and animations for the main examples are available for free 
as electronic supplementary material m- 

This paper is the outcome of the joint work of the Symbolic Computation group 
at RICAM (Linz): the problem of constructing linkages following a planar rational 
motion was proposed in our research seminar and after some time became its main 
topic. Each participant contributed to the final result according to his/her own 
background and skills (algebraic geometry, combinatorics, kinematics, etc). 

2. A FIRST EXAMPLE 

We want to illustrate the main ideas and techniques with an example. Consider 
the ellipse (x + l)'^ +4:y'^ = 1 in the plane, which admits the rational parametrization 

Our goal is to construct a planar linkage with rotational joints drawing the curve 
parametrized by ip that admits only one degree of freedom. This means that there 
is a specific link (on which we put the pen) that performs a motion along the ellipse 
as the linkage moves. By construction, the motion of this link is the composition of 
several rotations. We encode motions by univariate polynomials with coefficients 
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in the algebra K = ’C[rf\/{rf', irj + rji), where i denotes the imaginary unit. Elements 
of K are of the form z + r]w with z,w € C, and they are multiplied as follows: 

(z + 7)w) ■ (z'+ T/w') = (z z') + r/(zw'+ z'w). 

In our case, a translational motion along the ellipse is represented by the polynomial 

F(t) = (t^ + 1) + (it - 2). 

This means that the orbit of any point under this motion is a translate of the 
ellipse. To construct a linkage that realizes this motion, we want to factor P into 
linear polynomials, which correspond to revolute motions. However, one can check 
that this is not at all possible! On the other hand, as pointed out in Remark |4.8| 
we will see that for any R G M[t], the motion polynomial RP G K[f] describes the 
same motion as P. In the example, we can take R = + 1 so that RP admits a 

factorization into linear polynomials (see Theorem |5.15[ ): 

R{t) ■ P{t) = (t + i — r]i) ■ (t — i + ■ (t — i + ■ (t + i). 

Such a factorization is computed by the algorithm FactorMotionPolynomial (see 
the end of Section]^. It allows us to construct a linkage, in the form of an open 
chain (left part of Figure[^, whose links move according to the rotations represented 
by the linear factors (this is algorithm ConstructWeakLinkage at the beginning of 
Section]^. Since such a linkage has a high degree of freedom, we need to constrain 
its mobility by adding more links and joints (right part of Figure[^. This is achieved 
by an iteration of the so-called flip procedure, described in Sections and 




Figure 1 . An open chain and its extension to a linkage of mobility 
one, realizing the translational motion given by P{t). 


However, if we just want to draw the ellipse, we need not realize exactly the 
translational motion P{t): it is enough to find a motion for which the orbit of one 
point is the ellipse. As pointed out in Remark |4.9[ multiplying with a polynomial 
C G C[t] from the left does not change the orbit of the origin. In our case, we 
find that the polynomial CP with C{t) = t — i factors completely as follows (see 
Proposition 7.101: 


C(t) ■ P(t) = (t — i — \rii) ■ (t — 1+\r]i) ■ (t + 1 +Tji). 

This factorization gives rise to a slightly simpler construction, see Figure 
details concerning this example will be provided in Examples |4.12[ |5.16 


^ Further 
and |7.lH 
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Figure 2. The linkage that draws an ellipse. The same linkage 
is shown in different positions: t = 2 (white), t = ^ (light gray), 
t = 0 (dark gray), and t = —1 (black). 


3. Linkages 

In this section, we define a mathematical model for kinematic objects known as 
linkages. A linkage is a device constituted by rigid bodies, called links, connected 
by joints, which restrict the relative position of two neighboring links. In this paper 
we focus on planar linkages with revolute joints, namely linkages for which all links 
move in parallel planes and whose joints allow only rotations around a point. In our 
model we will not be concerned with the shape of the links — the only exception 
will be Section — and we will suppose that joints are the only constraints for 
the motion of the links. Thus we can represent a linkage by a graph whose vertices 
correspond to links, and where two vertices are connected by an edge if and only if 
the corresponding links are connected by a joint (for an example, see Figure]^. But 
the graph does not encode how the joints constrain the motion of the links, since it 
does not specify their (initial) position. So we have to add this extra information. 

Definition 3.1. A linkage with revolute joints is a connected undirected graph 
G = {V,E), together with a map p: E —K^, such that G does not have self¬ 
loops, i.e., all edges connect different vertices. The elements of V are called links, 
while the elements of E are called joints. We call two links neighboring if they are 
connected by a joint. For a joint e G E, the point p(e) is called the center of rotation 
of e. In the following, we will always assume that V is of the form {1,..., n} and 
that elements of E are given by unordered pairs {i,j} of elements i,j S V. 

Remark 3.2. An implementation of a linkage, where links are realized by line 
segments between the joints, also looks like a graph. Note that this graph is not 
the same as the one in Definition |3.H rather it is its dual. The lengths of the links 
are then given implicitly by distances between initial positions of joints. 
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Figure 3. A picture of a linkage (on the left) and the correspond¬ 
ing link graph (on the right). 


Next we explain how we describe configurations of a linkage, namely, the posi¬ 
tions of all its links at a given moment. There are essentially two ways to indicate 
the position of a link: we can give either its absolute position with respect to some 
frame of reference, or its relative position with respect to some other link. Both 
these notions will be encoded by means of isometries. 

Notation. We denote by SE 2 the group of direct isometries of i.e., maps that 
preserve distances and orientation (of the standard basis of M^): 

SE 2 := {cr: : cr is an isometry, deta = l}- 

For each link i there is a unique isometry ai G SE 2 taking it from its initial 
position to a current one; we say that ai is the absolute position of the link i. 
The relative position of a link i with respect to some other link j is defined to 
be (Tij = CTi o hence we have the relation Cij = ai^k o CTfc.i- If * and j are 
neighboring links, then (7^^ is a rotation around the point p = p(i,j), which is the 
initial position of the joint {i,j}- This follows directly, since o'i{p) = Cjip)- All this 
motivates the following definition of “virtual” relative positions; the name indicates 
that they do not take into account the constraints imposed by the linkage. 

Definition 3.3. Let L = {G,p) be a linkage, and let i,j G V be links connected 
by a joint e G E. The set of virtual relative positions VRP(i, j) of the link i with 
respect to the link j is the subgroup of SE 2 of rotations around the point p{e). 
Notice that VRP(i,j) = VRP(j,i). 

We define a configuration of a linkage L as a collection of virtual relative positions 
of links satisfying the following conditions: if (z, hi), (hi, h 2 ),..., (h^, z) is a directed 
cycle in G, namely a sequence of pairs of links connected by a joint, starting and 
ending at the same link, then the composition o ■■■ o — which gives 
the relative position of i with respect to itself — should be the identity. Loosely 
speaking, a configuration consists of the angles at all joints (counted twice), together 
with the (redundant) information p(e) for all e G E. 
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Definition 3.4. The set of configurations of a linkage L is defined to be 


Conf(L) := < {ak,i) G VRP(i, j) x VRP(j, i) : for every directed cycle 

^ {iJ}eE 

{i, hi), {hi,h 2 ), ..., {hs,i) in G, we have o • • • o ah„,i = idj 
where denotes the Cartesian product so that {ak,i) is a tuple of size 2\E\. 


This definition of configuration space is the main difference between our model 
and the ones appearing in the literature I31I1II31ITS1121]. 


Remark 3.5. The cycle conditions imposed in the definition of Conf(L) have 
the following consequence: suppose that i,j S V are two neighboring links; then 
(ijj), ij,i) is a directed cycle from i to i. Therefore, if S = {ck^i) is a configuration 
of L, the cycle condition imposes that aij o = id, implying that Cij = ckjl . 
Hence for every linkage L = {G,p), the projection 


H VRP(z,j) X VRP(j,f) 


n VRP(z,j) 

{i,j}eE 

i<j 


is a bijection when restricted to Conf(L) (and an isomorphism if we consider the 
projective structure on Conf(L) we will define soon), and similarly for every pro¬ 
jection that forgets exactly one among each pair {i,j) and (j, i). Still, we chose the 
above definition for Conf(L) because it does not fix an orientation of the edges, 
allowing us to deal with arbitrary directed paths (see Definition 3.71. 


Notice that for every pair of neighboring links i and j, the subgroup VRP(f,j) 
can be set-theoretically identified with the real projective line Pg. Under this identi¬ 
fication, every cycle condition imposed in Definition |3 .4| becomes a closed condition 
in the Zariski topology, since it is given by multihomogeneous polynomials. In this 
way Conf(L) acquires the structure of a projective subvariety of (P^) 


Definition 3.6. We define the mobility of a linkage L to be the dimension of the 
configuration space Conf(L) as a projective subvariety of (Pg) 

So far we only took into account the relative position of two neighboring links. 
For our purposes, namely to construct a linkage that follows a prescribed motion, 
we need to take one link (the “base”) as fixed, and consider the relative positions 
of all the other links with respect to the base. 


Definition 3.7. Let L = {G,p) be a linkage. Let E € Conf(L) and let i,j S U be 
links. Let (i, hi), (hi, /i 2 ),..., (hs,j) be a directed path in G from i to j — which 
exists, since by Definition |3 .1 1 the graph G is connected. Then we define the relative 
position of j with respect to i in the configuration E = (ufe,;) as 

RP(i,j, E) := CTj,III o ••• o e SE2. 

Notice that, because of the cycle condition, this definition is independent of the 
chosen path. We define the set of relative positions of j with respect to i to be 

RP(i,j) ;= {RP(i,j,E) : EeConf(L)} C SE2. 
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4. Motion polynomials 


As mentioned in Section our goal is to reduce the main problem to an algebraic 
one. We first introduce an algebraic setting for manipulating isometries, consisting 
in a noncommutative K-algebra K whose multiplication corresponds to the group 
operation in SE 2 . This is an instance of a general construction associating a Clifford 
algebra to each group of isometries SE„; see for example |26l Section 9.2]. Then we 
define the main object of our work, motion polynomials, as polynomials over K. 

Following m, one can embed SE 2 in the real projective space Pj| with coordi¬ 
nates xi,X 2 ,yi,y 2 , as the open subset 

(4.1) U := Pr\ {(a;i : X 2 : yi ■ y 2 ) : xl + xl= O}. 

Hence U is the complement of the projective line xi = X 2 = 0. The (right) action 
of an element {xi \ X 2 ■ yi ■ j/ 2 ) € Cl on a point (x, y) S is given by: 

xl - xl - 2 xiX 2 \ /Xiyi - X 2 y 2 \ 

2 xiX 2 xl-xl)\y) \xiy 2 +X 2 yi)/ 

Moreover, the group operation in SE 2 becomes a bilinear map: 

(xi : X 2 -.yi : 2 / 2 ) • (x'l : x '2 : y'l : y' 2 ) = - X 2 x '2 : Xix '2 + X 2 x'i : 

(4.3) xiy[ -f X22/2 + yi^'i - y2X2 ■ 

xiy2 - X2y[ + yix '2 + y2x'i) , 

where (xi : X 2 ■ yi ■ 1 / 2 ) and (x( ■. x '2 '■ y'l '■ y' 2 ) represent two direct isometries 
(T, cr' S SE 2 , respectively. 

For an easier handling of the multiplication in this embedding of SE 2 , we in¬ 
troduce the following notation: we write a representative (xi : X 2 : yi : 2 / 2 ) of an 
element of SE 2 as a pair (z, w) S C^, where z = Xi + 1 x 2 and w = 2/1 + * 2/2 (here i 
is the imaginary unit). Then Equation ( |4.3| ) can be rewritten concisely as 

(4.4) {z,w) ■ {z',w') = [z z',zw'+ z'w) 


(4.2) 


where z, z', w and w' are multiplied as complex numbers, and the bar (•) is complex 
conjugation. We can go further, by writing a pair {z,w) in the form z + rjw, and 
by postulating that rj satisfies the two relations: 

zr] = r]z for all z S C and rf = 0. 


Then the multiplication we obtain is exactly the one described in Equation (4.41: 
{z + yw) ■ {z'+ rjw') = z z'+ rj {z w'+ z'w). 


Definition 4.1. We define the M-algebra 

IK := C[r]]/{r]‘^,ir] + rji). 

Based on the previous discussions, we can identify elements of IK with elements 
of SE 2 . In this way the projective space Pj| in which we embed SE 2 can be thought 
as the projectivization P(K) of IK, considered as an K-vector space. 


Notice that the algebra IK is constructed in such a way that its multiplication is 
a lift of the group operation of SE 2 : this means that if two isometries cti, cr 2 G SE 2 
are represented by ki,k 2 G IK, then a 2 o cri is represented by fci • k 2 (remember we 
have a right action). More precisely, one can prove that K is isomorphic to the even 
subalgebra C^'''(0, 2,1) of the Clifford algebra (31(0, 2,1) of SE 2 ; see Section 9.2]. 
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Remark 4.2. Looking at the construction of the algebra K, we notice that the 
condition xf+x^ ^ 0 for points in Pj| representing isometries in SE 2 becomes z ^ 0 
when we consider elements of K. Moreover, the identity isometry is represented 
in Pg by the point (1 : 0 : 0 : 0), hence by any purely real element of K, namely 
elements of the form z + r]w with z G K and w = 0. From this we see that, given 
k = z+Tj w € K representing an isometry a G SE 2 , then k' = z — r]w represents 
since the product kk' equals jzp, which is purely real. 


The following lemma characterizes simple subgroups of SE 2 in a geometric fash¬ 
ion: it shows that both translational and revolute motions, i.e., families of transla¬ 
tions resp. rotations, correspond to lines in Pj|. 


Lemma 4.3. Let ^ C Pj| &e a projective line passing through the point (1 : 0 : 0 : 0), 


and define ijj = £ nU and X = £\U, where U is defined in Equation f.l Then: 


(i) if X has cardinality 1, then £u corresponds to a subgroup of SE 2 that 
consists of all translations along a fixed common direction; 

(ii) if X is empty, then ijj corresponds to a subgroup of SE 2 that consists of 
all rotations around a fixed common point. 


Proof. We analyze the two cases separately. 

Ad (i)- By hypothesis we have that (1 : 0 : 0 : 0) € f and (0 : 0 : a : &) S X C £ for 
some a, 6 S M, not both zero. Hence the line £ admits the parametrization: 


{X : 0 : apL : bp,), for {X : p) G 




Plugging this parametrization into Equation (4.21, we see that the elements of £u 
are translations by the vector u/x {a,b). 

Ad\(ii)\ Again by hypothesis we have that (1 : 0 : 0 : 0) £ moreover there exists 


a point in £ which is of the form {0 : a : b : c). Since X is assumed to be empty, we 
have a 0, and hence the line £ admits the parametrization: 


{X : ap : bp : cp), for [X : p) G Pj^- 

We use Equation ( |4.2| to compute the fixed points of an arbitrary element of £. 
A direct calculation reveals that we always get the point {—‘^f 2 a,bj 2 a), which is 
independent of A and p. Thus every point of £u represents a rotation around it. □ 


The description of revolutions around a point given by Lemma |4.3| enables us to 
compute effectively the configuration curve of a linkage, as shown in Example |4.4[ 

Example 4.4. Let us consider the linkage L whose graph is depicted in Figure]^ 
The linkage L is given by 4 links and 4 joints, and the map p is determined by: 


P(l,3) 

= Ml 

= (0,-3/2), 

P(3,4) 

= M3 

= (-VaVs), 

P(2,4) 

= U4 

= (-1,13/2), 

P(l,2) 

= M2 

= (-VaV2). 


We compute the configuration set of L and its dimension as a projective variety. 
Since we have 4 joints, it will be a projective subvariety of (Pr)"^ (here, rather than 


Conf(L), we are considering one of its projections, but as noticed in Remark 3.5 


we get an isomorphic object). Lemma 4.3 ensures that the 4 subgroups of SE 2 of 
rotations around the points Ui — which constitute the subgroups of virtual positions 
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Ui U4 

Figure 4. Graph of a linkage constituted by 4 links and 4 joints, 
a so-called closed 4ii-linkage. The edges are labeled by points Ui 
in that are the images of the edges under the map p. 


between two neighboring links as in Definition 
with parametrizations: 


3.3 


correspond to 4 lines C Pj| 


center u, 


parametrization of line (.i 


( 0 ,- 3 / 2 ) 

(-3/4, 9/2) 

(- 1 , 13 / 2 ) 
(-1/4, 1/2) 


(Ai : pi : —3pi : 0) 

(A 3 : 2p3 ■ 18/43 : Sps) 
(A4 : Pa '■ 13/44 : 2/44) 
(A2 : 2/42 : 2/42 : P2) 


where (Ai : pi) G Pg for all i € {1,...,4}. Recall from Definition 3.4 that, in 
this case, a configuration for L is a 4-tuple (o'i_ 2 ) o' 2 , 3 , (T 34 , cr 4 ^i) of direct isome¬ 
tries that satisfies the cycle condition cri _2 o <72,3 o < 73,4 o < 74^1 = id. Each of the 
isometries ak^i gives a point on one projective line £i, and in our projective model 
of SE 2 composition corresponds to multiplication according to Equation (4.31. Un¬ 


der these identifications, the composition cri ^2 ° <^ 2,3 o ( 73^4 o (74 corresponds to 
a point {Fi : F2 : F3 : Fa) G Pj|, where all Fi are polynomials in the variables 
(Ai : pi), ..., (A 4 : Pa)- Noticing that the identity element of SE 2 is represented by 
the point (1 : 0 : 0 : 0) € Pj|, one realizes that the cycle condition is equivalent to 


rk 


Fi F2 F3 F4 
10 0 0 


= 1 . 


Hence, as a subvariety of (Pr)^, the configuration set Conf(L) is the zero set of 
the polynomials F 2 ,F 3 ,Fa- A computer algebra computation shows that this is a 
one-dimensional variety with two components. 


We introduce now one of the main concepts of this paper, namely motion polyno¬ 
mials. Intuitively, a motion can be described as a curve in the space of direct isome¬ 
tries. In our case, we want such a curve to be defined by a rational parametrization. 

Definition 4 . 5 . Let Xi, X2,Yi,Y2 G K[t] be polynomials such that A/ -I-A| is 
not identically zero; denote by U(Ai, A 2 , Yi, I 2 ) the set of their common zeros 
in K. The map (/: IR\ U(Ai, A 2 , Yi, Y 2 ) —^ Pr defined by Ai,A 2 , Yi,Y 2 is called 
a rational motion. We will use the notation (/: K ---» Pr to mean that (f> is not 
defined everywhere on M. Thinking of SE2 as an open subset of Pr, the condition 
on A4 -I- A| ensures that for all but finitely many t gM. we have (j){t) G SE2. 

Definition 4 . 6 . Let be a rational motion given by (Ai, A2, Yi, Y2). The polyno¬ 
mial P{t) = Z(t) -G rj W(t) G ]K[/], where Z = Xi-G 1X2 and lY = Yi -|- 4 Y 2 , is called 
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a motion polynomial encoding the motion (j). The polynomials Z and W are respec¬ 
tively called the primal and secondary part of P, denoted by pp(P) resp. sp(P). 


Next we want to connect motion polynomials with rational plane curves. For this 
purpose it is helpful to rephrase the action of SE 2 on shown in Equation (4.21 in 
the following way: we identify a point {x, y) G with the element u = x + ly G C; 
if a G SE 2 is represented by an element z + rj w G then cr sends u to 

(4.5) 


uz 


■ zw 


Using this formulation of the (right) action, it is easy to prove the following result. 


Proposition 4.7. Let ip: K —> be a rational parametrization of a real curve, 

which means that ip is of the form 



for some real polynomials /, g and h. Then the orbit of the origin under the motion 
given by the motion polynomial P = h + g (f + ig) is exactly the image of p. 


Let Pi and P 2 be two motion polynomials. Because of the algebraic properties 
of K, the polynomial P — P 1 P 2 defines a motion that is the composition of the mo¬ 
tions determined by P 2 and Pi. More precisely, for every t G M., the isometry P{t) 
is the composition of the isometries P 2 (t) and Pi{t). This shows how important 
the factorization of a motion polynomial is in our framework: it provides a decom¬ 
position of a motion into simpler ones. 

We saw in Lemma 4.3 that rational motions whose image is a line in Pj| are 
translational motions or revolutions. Hence linear motion polynomials encode this 
kind of motions. It follows that, if we are able to factor a motion polynomial into 
linear ones, then we can decompose a rational motion into revolutions and trans¬ 
lational motions. Later we will consider a suitable subclass of motion polynomials, 
so that only the first situation occurs. 


Remark 4.8. Let (f)'. M Pj| be a rational motion and let P{t) G K[t] be the 
corresponding motion polynomial, so P = Z + g W with Z,W S C[t]. Notice that, 
although Z and W have coefficients in C, the polynomial P determines a real curve 
in n. namely a curve described by real equations. Hence, if P G K[t] is a nonzero 
real polynomial, one can check that PP G K[t] provides the same curve in Pj| as P. 

Remark 4.9. Let again (f): K — + Pjg be a rational motion and let P{t) G IK[t] be 
the corresponding motion polynomial. For any C G C[f] the orbits of the origin 
under the motions given by P and CP are the same since C represents a revolution 
around the origin — recall that CP acts on by first applying C and then P. 


Example 4.10. The rational motion given by the polynomial f -I- ?7 is a vertical 
translational motion. Indeed, its image in Pj| is parametrized by (A : 0 : /r : 0), so 
that we are in case (i) of Lemma 4.3 In contrast, the motion polynomial P(t) = 
t + i gives a revolution around the origin (0, 0). This is an instance of case (ii) of 
Lemma 4.3 using the parametrization (A : /i : 0 : 0). 

Example 4.11. Consider the following product of two linear motion polynomials: 


(t-l-i) • {t — i + g) = {t'^ + 1) + g {t — i). 

' -V-^ ^-^ '-V-' 

rotation around rotation in the opposite this is a translation, since 

a point sense, around another point there is no imaginary part in t^ + 1 
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For a fixed t we get a translation by the vector —1). As t changes, this vector 

describes a circle, hence we get a translational motion along a circle (see Figure]^. 



Figure 5. The trace of the point (0,0) under the circular transla¬ 
tion given by the motion polynomial + 1 ) + rj (t — i) is depicted 
via a dashed line. The action of the motion on the black triangle 
highlights the fact that we have a purely translational motion. 


In Example |4. 11 1 we saw that the factorization of a motion polynomial into linear 
polynomials provides a decomposition of the described motion into revolutions. 
However, Example |4.12| shows that this is not always possible. 


Example 4.12. Let us consider an elliptic translational motion: 

P{t) = {t^ + \) + Tj {at — hi), a,bGR. 


We try to factor P into two linear polynomials Pi, P 2 & K[t]. First notice that the 

primal parts of Pi and P 2 should be factors of -f 1 ), so they have to be of the 

form t Pi. Thus we have only two possibilities: 

f PP(-Pi) = t + i I pp(Pi) = t-i 

1 PP(-P 2 ) = t-i \ pp(P 2 ) = t + i 

By a direct computation one can prove that none of those two choices gives a 
factorization if a ^ 6 , thus P cannot be factored into linear polynomials. 

On the other hand, we show now that if we multiply P by a real polynomial R 
we can achieve a factorization. Recall that, by Remark |4.8| the polynomials P 
and RP describe the same motion. In this case we take R — t^ + 1. Hence we need 
four linear polynomials Pi,..., P 4 to factorize RP, and again pp(Pi) = t ± 1 . We 
make the following ansatz: 


(4.6) 


Pi 

P2 


Imposing Pi • • • P 4 


= {t-i)+r]Wi, P3 = {t + i)+r]W3, 

= {t + i)+r]W 2 , P4 = {t — i)+r]W4. 

= RP gives a linear system in the Wi, namely: 


(4.7) 


( wi + W 2 + W 3 + W 4 = a, 

\ Wi + W 2 — W 3 — W 4 = b. 


Each solution of the system (4.7) gives rise to a factorization of the polynomial RP. 
If, instead, we take the following ansatz: 


Pi = {t + l)+f]Wi, P 3 = {t + l)+f]W3, 

P2 = {t-i)+r]W2, P4 = {t-i)+r]W4. 


one can check that the corresponding linear system does not admit any solution. 
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In the next section, we consider the factorization problem in a more systematic 
way. We will see that, once we restrict to a certain class of motion polynomials, all 
situations can be treated as we did in Example |4.12[ in general a motion polyno¬ 
mial P cannot be factored into linear polynomials in K[t], but it is possible to find 
a real polynomial R such that RP can be factored. 


We end this section by defining precisely what we mean by saying that a linkage 
realizes a given rational motion. Recall from Definition |3.7| that, given two links i 
and j of a linkage L, we can define the set RP(z, j) of relative positions of the link j 
with respect to the link i, and this is a subset of SE 2 . Furthermore, recall that the 
set Conf(L) of configurations of L has the structure of a real projective variety. We 
are going to use some concepts and results coming from real algebraic geometry: 


(i) every projective space is a real affine variety, namely there exists an 


(ii) 


(iii) 


embedding 




for some k G N, and C(lPffi) is the zero set of a 


collection of real polynomials; see for example |3J Theorem 3.4.4]; 
a subset X C is called semialgebraic if it can be described as the set 
of points satisfying a disjunction of conjunctions of polynomial equalities 
and inequalities; see for example O Definition 2.1.4]; 
a subset X C Pg is called semialgebraic if its image under any embedding 
^: Pg as in 


(i) is a semialgebraic set in the sense of (ii) 


Lemma 4.13. Let L be a linkage and let i,j be links of L. Then, under the 
projective embedding 0 /SE 2 , the set RP(z, j) becomes a semialgebraic subset 0 /SE 2 . 


Proof. Let (z, hi ),..., {hs,j) be any directed path from z to j. Consider the map 
F: Conf(L) —SE 2 sending a configuration E = {o'k,i) to ' 'Ocrh„,j- One can 

check that RP(z,j) coincides with F(Conf(L)). Moreover, if we write the map F 
in terms of the projective coordinates of Conf(L) and SE 2 , then F is given by real 
polynomials. Hence, because of a general result in real algebraic geometry (see for 
example |3J Theorem 1.4.2]), the set RP(z,_)) is semialgebraic. □ 


In the following definition, and later in this paper, the word “component” will 
always stand for “irreducible component” with respect to the Zariski topology. 

Definition 4.14. Let L be a linkage and let K Pj| be a rational motion. 
We say that 

(i) L weakly realizes the motion (j) if there exist links z and j of L such that 
image(())) C RP(z, j). 

(ii) L strongly realizes the motion (j) if there exist links z and j of L such that 
image(0) C RP(z,j) and image((/)) is a component of RP(z,j), where (•) 
denotes the Zariski closure in RP(z, j). 


Remark 4.15. Notice that in case (ii) of Definition 4.14 we allow the given variety 
of relative positions to have several components, and we ask that only one of them 
coincides with the Zariski closure of the image of the given rational motion, which is 
irreducible by construction. This situation occurs in the example given on the right 
of Figure]^ since the variety of relative positions between base and end effector has 
two components (see Lemma 6.6). 
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Figure 6 . Two linkages that both realize the motion give n by 
(t + i — 377) • (t + 2z + 77(18 + Sj)) , discussed in Example |4.4l The 


depicted position is reached at t = —2. The linkage on the left has 
two degrees of freedom and thus weakly realizes the given motion 
as the relative position between base (in black) and end effector 
(in white). The linkage on the right has mobility one and provides 
a strong realization of this motion. 


5. Factorization of motion polynomials 

Motivated by the considerations from the previous section, we now address 
the problem of factorizing a motion polynomial P = Z + rjW G IK[t] into lin¬ 
ear factors. The main results are Theorem |5.15| and the corresponding algorithm 
FactorMotionPolynomial. To achieve them, we have to go through some techni¬ 
calities, which however will not be needed in the rest of the paper. 

We restrict our attention to monic polynomials, i.e., to polynomials P whose 
leading coefficient is 1, which implies that degW < degZ. Our goal is to write a 
monic P S IK[f] in the form P = Pi • • • P„, where n = deg P and deg Pi = 1 for all i. 
If P = Z + r]W and Pi = Zi + r]Wi, then Z = Z\ - ■ ■ Zn- Hence the primal part of 
each Pi has to be one of the linear factors of Z G €.[t]. We start by a few definitions 
and results characterizing monic motion polynomials that split into linear factors. 

Definition 5.1. Let z = {zi ,..., Zn) G C". We define the following polynomials: 

Qi{z) := (t - W) • • • (f - WZi)it - Zi+i) Zn) for 7 e {1,..., n}. 

Hence we have that degQi(z) = n — 1 for all i. 

Remark 5.2. Let Z G C[f] and let z = {zi, ..., z„) be a fixed permutation of the 
roots of Z. Then the polynomials Qi{z) have a non-trivial gcd if and only if the 
polynomial Z admits a pair of complex-conjugate roots a and a (this includes the 
case of multiple real roots). 

Lemma 5.3. Let P = Z + i^W G K[f] be a monic motion polynomial, and let 
z = {zi,..., Zn) be a fixed permutation of the roots of Z over C. Then P admits a 
factorization P = Pi • • • P„, where Pi{t) = (t — Zi) + rjWi with Wi G C, if and only 
if W lies in (^Qi{z), ..., (5n(z))(;, i-e., in the C-linear span of Qi{z), ..., (5„(z). 
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Proof. Suppose that P G K[t] admits a factorization into linear factors: 

n n ^ \ / ^ \ 

P{t) Zi + VWi) = Y[{t-Zi){t-Zj)jwk. 

i—1 i—1 k—1 ^ j—1 ^ ^ ^ 

" -v-" 

= Qfc(z) 

Such a factorization exists if and only if we can choose wi,..., Wn G C such that 
J2k=i'^kQk{z) matches the prescribed 7?-part of P{t) = Z{t) + riW{t), namely if 
and only if Vh e (( 5 i(z),... ,Q„(z))j,. □ 

Lemma 5.4. Let P G K[t] be a monic motion polynomial and suppose that Z = 
pp(P) has no pair of eomplex-eonjugate roots. Then for every permutation z = 
(zi,..., Zn) of the roots of Z, the polynomial P admits a faetorization P = Pi ■■■ Pn 
with pp(-Pi) = t — Zi for all i. 


Proof. We fix a permutation z of the roots of Z, and we denote Qk = Qfe(z) for 
all k. From Lemma [5.3| we know that P = Z + rjW admits a factorization if and 
only \iW G {Qi, ..., Qn)^- Clearly, this is always possible (for arbitrary W) if the 
following matrix £ (^nxn jg non-singular: 


Mn 


(t')Qi 


{t^)Qn 


V(r-i)Qi ••• {t'^-^)Qn) 


Here {P)Qk denotes the coefficient of T in the polynomial Qfc- Notice that the 
matrix entries are, up to sign, elementary symmetric polynomials in the Zi and zi. 
We now exhibit that the determinant of Mn is nonzero under the assumptions on 
the roots of Z. Indeed, we find that 


(5.1) det(M„) = (zi-Zj), 

which follows, as a special case, from Lemma 3 in jTS] (by substituting Ak = —Zk, 
Bk = —Zk-i, and by extracting the coefficient of ^ similar deter¬ 

minant evaluation is given in m where the Zi appear without conjugation. Since 
det(M„) is very much reminiscent of the Vandermonde determinant, it is no surprise 

□ 


that (5.1 1 can also be proved in an analogous fashion. 

,Zn) G C" such that gcd(^Qi{z),... ,Qniz)) = 1. 


Corollary 5.5. 

Then 


Let z = (zi,..., 


(Qi(z), ■ • ■ ,Qn(z))c = {W G<C[i\ :deg(W)<n}. 


Proof. The gcd condition implies that z has no pair of complex-conjugate roots, as 
noticed in Remark 15.21 Then the claim follows from Lemma [5)3l and Lemma□ 


Since the gcd of the polynomials Qi{z) introduced in Definition 5.1 plays an 
important role for the factorization of motion polynomials, we make some effort to 
describe it precisely (Lemma 5.6 Corollary |5.7[ Proposition 5.10); for this purpose 
the combinatorial notion of matching will be introduced. At the same time we 
aim at a nicer characterization of those secondary parts in Lemma |5.3| that en¬ 


sure the existence of a factorization (Proposition 5.12). To simplify the proofs of 
Propositions |5.10| and |5.12[ we formulate the somewhat technical Lemma |5.9| 
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Lemma 5.6. Let a, /3 € C with a ^ f3 ^ a, and let z, = (zi,..., Zn) G C" and z = 
(zi,..., Zj-i, /3, Zj,..., Zn) G for some j G {1,..., n + 1}. Then the multiplic¬ 

ities of the roots a and a in gcd(Qi(z), ■ ■ •, Qn(z)) and in gcd((5i(z),..., Qn+i(z)J 
are the same. 

Proof. First observe that Qfc(z) and Qe{z) have the same multiplicities of a and a, 
where t = kiik < j and £ = fc+lifA:>j. On the other hand, the polynomial Qj(z) 
has no such counterpart among the Qkiz), but its multiplicities of a and a are 
greater than or equal to those in its (one or two) neighbors Qj-i(z) and Qj^i(z). 
The claim follows. □ 


Corollary 5.7. Let x G C™ and y G C" with Xi ^ Vj ^ xi for all 1 < i < m 
and 1 < j < n. If z G C™"*"" is an arbitrary interlacement of x and y, i.e., 
there are l<ii<---<im^^ + n and 1 < ji < ■ ■ ■ < jn < m n with 
{zi,... n {ji,... ,j„} = 0 such that Zi^ = Xk, 1 < k < m, and zj^ = yk, 
1 < k < n, then 

gcd(Qi(z),..., Qm+n{'Z‘)) = gcd(Qi(x),... ,Q„(x)) • gcd(Qi(y),... ,Q„(y)). 
Proof. This is a direct consequence of Lemma |5.6| □ 


Definition 5.8. Let z = (zi,..., z„) G C”. A set 

Af C {{i,j):l<i<j<n and Zi = 

is called a matching of z if for all {ii,ji), (* 202 ) G M we have i\ ^ ii and ji ^ j 2 - 
(Note that for each z G C" there is a bipartite graph with at most 2n — 2 vertices 
such that the matchings in that graph and the matchings of z are in bijection.) 


Lemma 5.9. Let z = (zi,..., z„) G C", let G = gcd(Qi(z),..., (3„(z)), and let 
g = degG. For a matching M of z define H{M) := Y[{i j)eMi^ ~ ^ 3 )- Then the 
following statements hold: 

(i) For any matching M of z we have H{M) \ G. 

(ii) There is a matching M of z such that H{M) = G. 

(iii) There are indices 1 < zi < • • • < *n-g < n such that, setting z = 
(zq,... ,Zi„_g), we have Qk{z) = Qi,.{z)/G for all 1 < k < n — g and 
gcd(Qi(z),...,Q„_g(z)) =1. 

Proof. For|(i)|we show that H{M) divides Qk{z) for any matching M of z and for 
any k G {1,..., n}. Indeed, by rewriting 



j>k j<k 


the claim follows, since Qki^) has precisely the roots ., Zfe_i, Zk+i, ■ • ■, Zn. 

For (ii) and |(iii) note that the statements trivially hold for G = 1, so we assume 
now that G 7 ^ 1. Using Corollary |5.7| the proof is reduced to the case where 
z G {a. O'}” for some a G C. If a G M then we have Qi(z) = • • • = Qn(z) = 
{t — a)"“^ = G and we can choose M := {(1,2), (2, 3),..., (n — 1, n)} and ik = k 
for 1 < fc < n. Hence from now on we assume a a. 

Both assertions |(ii)| and (iii) are proved in parallel by an inductive argument 
on the degree of G, so in the following we show how to reduce the given scenario 
to one with smaller g. The main difficulty here is to pick an index j such that 
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deleting Zj from z yields a gcd of smaller degree. We distinguish three cases and 
assume w.l.o.g. that = a (case A) resp. Zi= a (cases B and C). 

Case A: There exists an index j such that ^^{zj,..., Zn) < ^a{zj, ■ ■ ■, Zn), 
where denotes the number of occurrences of a. In this situation we set z = 
{zi,..., Zn-i) and G = gcd((5i(z),... ,Qn-i(z)). By the induction hypothesis 
there exists a matching M for z such that H{M) = G. Because of the assumption 
on the number of a’s and a’s, there must be an index i € {j,... ,n — l} with Zi =a 
that does not occur as the first entry of any pair in M. Hence M := M U {{i, n)} 
is a matching for z and we have H{M) = {t — a)H{M). Since by (i) we have 
H[M) I G, we obtain {t — a)G \ G. Next observe that Qk{i) = Qk{'^)/{t — Zn) for 
all 1 < fc < n and hence 

G = gcd((3i(z),..., Q„_i(z)) = gcd((5i(z),..., Q„_i(z)). 


From the definition of G it follows that G | {t—a)G, and thus G = {t—a)G = H[M). 

Case B: There exists an index j such ^^izi, ■ ■ ■, Zj) < ^a{zi, ■ ■ ■, Zj). The 
reasoning in this case is completely analogous to case A, setting z = {z 2 , • ■ •, Zn). 

Case C: Neither case A nor case B applies, which means that for all indices j we 
have ^a{zi^ ■ • ■ 5 Zj^ > • • • j Zj) and • • • ■> Zn) ^ 4fa(^Zjj ■. •, Zn^. Note 

that in this situation we automatically have zi = Zn, which w.l.o.g. we assumed to 
be a. Now there exists no entry in z whose removal would yield a gcd of degree 
deg G — 1. We circumvent this problem by setting z = (zi,..., Zj_i, Zj+ 2 , ■ • ■, Zn) 
where j is the largest index such that zj = a; in particular j < n. The rest 
is analogous to the previous reasoning. As before, we may assume that there 
exists a matching M of z with H{M) = G = gcd(Qi(z),..., ( 3 n- 2 (z)). Since 
^a{zi,... ,Zj) > #q(zi, ... ,Zj) there exists an index i < j with Zi = a such that 
i does not appear as the first entry of any pair in M. In order to extend M to 
a matching for z we need some relabeling; for this purpose we define r(a) = a if 
a < j and T(a) = a + 2 otherwise. Then 

M := {(T(a),r( 6 )) | (a,5) e M} U {(*,j),(j,j + 1)} 

is a matching for z and H[M) = {t — a){t — a)H{M); hence (t — a)(t — a)G \ G. 
A simple calculation shows that Qk(z) = Qk{z)/(^{t — a){t — a)) for k < j and that 
Qk{z) = ( 5 fc+ 2 (z)/((t - a){t- a)) for k > j. Thus we get 

^ = (^t - a){t - a) ■ • ■ ’ Qj+ 2 (z), ■ ■ ■, Qn(z)) • 

Again, this implies G = {t — a){t — a)G = H{M), which concludes the proof. □ 


Proposition 5.10. Let z = (zi,...,z„) G C" and let M be a matching of z, of 
maximal size, and let Qi{z ),..., (5n(z) be as in Definition 5.1 Then we have 

G :=gcd(Qi(z),...,Q„(z)) = {t - Zj) 

(*j)GM 

(where the gcd is assumed to be a monic polynomial). 


Proof. In Lemma 5.9 part |(i)| it was shown that the product on the right-hand side, 
which we denoted by H{M), divides G for any matching M. For the other direction 
we have to argue that the maximality of M implies that G | H{M), or equivalently, 







18 


GALLET, KOUTSCHAN, LI, REGENSBURGER, SGHICHO, AND VILLAMIZAR 


that G \ H{M) implies that M is not maximal. From Lemma 5.9 part (ii) it 
follows that there exists a matching M' of z such that H{M') = G. Note that 
\M'\ = deg(G'). Thus, if G -j" H(M), which means that deg{H(M)) < deg(G), then 
M is not maximal because |M| < \M'\. □ 


Example 5.11. Let Z = (t — aY{t — aY'^^ for some a S C \ K. In the following 
table we consider different permutations z G of the roots of Z; to each 

permutation we give the gcd G = gcd{Qi{z),... ,Qn{z)) and a matching M of 
maximal size that witnesses G: 


permutation G M 

{a,...,a,a,...,a) {t-aY {(1, r + 1 ), (2, r + 2),..., (r, 2r)} 

{a,...,a,a,...,a) {t-aY {(1, ?'+ 2), (2, r + 3), ..., (r, 2r + 1)} 

{a,a,a,a,... ,a,a) {t — aY{t — aY {(1, 2), (2, 3),... , (2r, 2r + 1)} 


The cases displayed above are the extreme ones: indeed, it is easy to see that 
r < deg(G) < 2r. Moreover, for any G = {t — ay{t — aY with 0 < i,_;/ < r 
and i + j > r there exists a permutation z that produces exactly this gcd G. 
These considerations lead to interesting combinatorial questions — e.g., how many 
permutations of a given set of roots are there that produce a prescribed gcd — 
which, however, are irrelevant for our purposes. 


Proposition 5.12. Let z = ( 21 ,..., z„) G C" and W G C[t] with deg IT < n. 
Then 


IT G (Qi(z),... ,Q„(z))j, IF G (Qi(z),...,Q„(z)) •C[f], 


where (Qi(z),..., (5„(z)) • C[t] is the ideal ofC[t] generated by (5i(z),..., Q„(z). 


Proof. The claim is equivalent to the following statement: 

(5.2) (Qi(z),...,Q„(z))j, = (Qi(z),...,Q„(z)) •C[t] n C[t]<„, 

where C[f]<„ denotes the set of complex polynomials of degree less than n. Note 
that the contain ment “C” is trivial. Let G = gcd{Qi{z),... ,Qn{z)). If G = 1, 
then by Corollary 5.5 we have (Qi(z),..., Qn{z)')^ = C[t]<„, which implies Equa¬ 


tion (5.21. On the other hand, if G is non-trivial. Equation (5.2) is equivalent to 
(5.3) 


{Qi,...,Qn)c = (Qi,...,Q„)-Cit] n C[t]<„_ 3 , 


5.9 


part (iii) 


there exists z = 


where Qi = Qi{z)jG and g = degG. By Lemma 
( 2 jj, .. such that gcd((3i(z),... ,Q„_g(z)) = 1, and (5fe(z) = for all k. 

Hence, by Corollary 5.5 we have (Qi(z),.. ■, Qn-g(z})^ = C[t]<„_g. This implies 


Equation ([^, since (Qi(z),... ,Q„_g(z))p C (Q^,..., Q„)^. 


□ 


Recall that our goal is to decompose a given motion into a sequence of revolu¬ 
tions, by factorizing a corresponding motion polynomial. For this purpose we have 
to restrict the domain of motion polynomials we are working with: for example, 
there is no hope to write a motion whose orbits are unbounded as the composition 
of revolutions, since the orbits of revolutions are always bounded. In particular, if 
we consider a rational bounded curve, then it can be given by a parametrization 
ip = {f /h,g/h) for some real polynomials /, g and h with deg h > maxjdeg /, deg g}. 
The boundedness of the curve implies that h has no real roots. The limit of (f{t) 
for t —>■ (X) gives a point in K^, which can be translated to the origin. Then we even 
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have degh > max{deg/, degg}. If we further assume that h is monic, then also 
the motion polynomial P = h + rj {f + ig) is monic, and by Proposition 4.7 it traces 
the image of ip. This motivates the following definition. 


Definition 5.13. Let P = Z + rjW he a motion polynomial in K[t]. We say that P 
is bounded if it is monic and if Z does not have real roots. Notice that the motion (j) 
given by a bounded motion polynomial is a map defined for all t € K; hence for 
such a motion we can use the notation (j): K —Pj| instead of K Pj|. 

Remark 5.14. A bounded linear motion polynomial P represents a revolution 
around a point: if P{t) = t — k, then t — k € K. represents a rotation around a 
point Q that depends only on k. This follows from Lemma [4.3[ 


We are now ready to state and prove the main result of this section: 


Theorem 5.15. Let P S IK[t] be a bounded motion polynomial. Then there exists 
a real polynomial R S M[t] such that RP can be factored into linear polynomials. 

If P — Z + rfW and R = gcd{Z, Z), then the smallest such real polynomial is 
R = R/ gcd{R, WW). 

Proof. If P is of the form P = SP with S G K[t]\K and P G K[t], then we can apply 
the theorem to P obtaining R G K[t] such that RP factors; thus also RSP = RP 
factors. Hence we may assume that P does not have non-constant real factors. 

Using Lemma |5.3| and Proposition |5.12| the proof is reduced to showing that 
there exists a permutation of the n roots of RZ (counted with multiplicities) such 
that the corresponding G = gcd((5i,..., Qn) divides RW. To this end write 


Z = r^>Si> 0, > 0, 

where € C are pairwise conjugate-distinct, i.e., ^ and Ui of for i j. 

Note that in the special case Si = • • • = = 0 we get i? = 1, in agreement with 

the fact that in this case already P itself factors by LemmaOn the other hand, 
we have always deg(i?) < deg(P), and equality is attained when ri = Si for all i 
and gcd{Z, W) = 1. 

Define Ri := {t — ai){t — of) and write R = JliLi ■ Next let Ui and vi be the 
multiplicities of at and of in W, i.e., 

h 

W = W-l[(^{t-a,nt-cr,r^ 

i=l 


for some polynomial W G C[t] such that gcd(IU, i?i • • ■ Rh) = 1- By introducing 
the quantity mi := minjsi, Ui -I- Uj}, the polynomial R can be written as 


R = 


R 

gcd{R, WW) 


n«r- 
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and hence 


RZ = - oti 

Z=1 




RW = W 

To construct an admissible permutation, the roots are arranged as follows: 


Cli, ■ ■ ■ 5 • 1 CI7, OL 




. . . , 0^2 ^ . 


Si—min{si,i!i} Ti'\-Si — 'mi Si —min{si,Ui} 

The assumption on P ensures that if Si > 0 then Ui = Q ov Vi = 0; indeed, if Si,Ui, Vi 
were all strictly positive, then i?j would be a common real factor of Z and W. 
Then = min{si,Ui} + min{sj,Ui}, and hence the previous multiplicities agree 
with those in RZ. 

The final permutation is obtained as an arbitrary interlacement of all arrange¬ 
ments for 1 < z < /i. Proposition |5. 10 then implies that 

h 

Since — min{si, Ui} < Ui + Si — rrii and — min{sj, < Vi + Si — rrii, we conclude 
that G divides RW. □ 


Algorithm 1 FactorMotionPolynomial 

Input: P = Z + r]W G K[t] a bounded motion polynomial such that Z and W 
have no common factor in R[f] \ K. 

Output: a polynomial R G M[t] and a tuple (fci,..., fc„) of elements of K such that 
(t - ki) ■ ■ ■ (t - kn) = R{t) ■ P{t). 


1 : Factor Z{t) over C, obtaining Z = ~ ctiY'it — of)®’ with Vi > Si> 0. 

2 : Initialize q — empty (the empty tuple). 

3: For i = 1,... ,h Do 

4: Set u-i — muxj ^i')^ I uud Uj — muxj | 

5: Set rrii = min{si, Ui + Vi\ 

6 : Set U) — ^ G-i , . . . , Gi , Gi , ■ ■ ■ ; Gi^ Gi , . . . , Gi ^ . 

Sj —min{si,i!i} ri-\-Si—mi —min{si ,Wi} 

7: Set q = concatenate(q, w). 

8 : End For 

9 : Set i? = n!Ll((^ - «i)(^ - 

10 : Set n = length(g). 

11 : Set Qj = nCi nr=j+i(^ - aU j e {1,... ,n}. 

12 : Compute {wj}^^i such that RW = using linear algebra. 

13: Set kj = qj — qwj for all j G {1,..., n}. 

14: Return (i?, (fci,..., A:„)). 


The proof of Theorem |5.15| immediately gives rise to a factorization algorithm 
for motion polynomials, called FactorMotionPolynomial. It produces one single 
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factorization for a given polynomial P, although there could be many for the follow¬ 
ing reasons: first, there may be several admissible permutations of the roots of RZ, 
and second, the choice of the coefficients {wj} in Step 12 may not be unique, since 
the polynomials {Qj} need not be C-linearly independent. 


Example 5.16. We can review Example |4.12| in the light of Theorem |5.15| and 
the algorithm FactorMotionPolynomial. We consider a particular instance of the 
collection of motion polynomials 

P{t) = + 1) + r] {at - bi), 


namely we choose a = i and b = 2t, obtaining the motion we presented in Section]^ 
If we run FactorMotionPolynomial with P as input, we obtain R = + 1, and if 

we fix the permutation {i, —i, —i,i) of the roots of pp(i?P), as in ( |4.6[ ), we obtain 
a two-dimensional family of factorizations into linear factors, i.e., the solutions of 
the linear system (4.7), out of which we choose the factorization used in Section]^ 


6. The flip procedure 


The procedure described in this section is crucial for the construction of linkages 
with mobility one. It is inspired by similar techniques — involving the interchange 
of factors in a factorization of a quadratic motion polynomial — used in [3] and m- 

Lemma 6.1. Let ki,k 2 S IK &e such that pp(fci) ^ pp(fc 2 ). Then there exists a 
unique pair (^ 3 ,^ 4 ) G such that: 

(i) pp(fc 3 ) = pp(fc 2 ) and pp(fc 4 ) = pp(fci); 

(ii) {t — ki){t — k 2 ) = {t — k^){t — ki) as polynomials in IK[t]. 


Proof. Let us suppose that ki = Zi + rjWi for j G {1,2}. We make the ansatz 
ks = Z 2 + r]W 3 and ^4 = zi + pw 4 , where W 3 and W 4 are elements of C to be 
determined. Then condition (ii) is equivalent to the linear system 


/ 1 1 \ / W3\ ^ f Wi +W2 \ 

V-^l Z2J \W4j \ZiW2+Z2WiJ 


which has a unique solution, since by hypothesis the determinant of the matrix is 
different from zero. □ 


Definition 6.2. Let ki,k 2 G K be such that pp(fci) ^ pp(fc 2 ). Then we define 
Flip(A:i, ^ 2 ) •= (ks, ^ 4 ), where (^ 3 , k 4 ) G is the pair from Lemma 6.1 


Definition 6.3. Let k = z + rjwGK with z G C \ K. Then k represents a rotation 


around a point that we denote by x(/c). Using Equation (4.5) one obtains that 

X(k) = 

z — z 


Definition 6.4. Let fci, ^2 G K be such that pp(fci) ^ pp(fc 2 ) and pp(fci), pp(fc 2 ) G 
C \ K. We define the flip linkage of ^ 1,^2 via the algorithm FlipLinkage below. 
The graph of this flip linkage is shown in Figure]^ 

Lemma 6.5. Let ki,k 2 G K with pp(fci) ^ pp(fc 2 )- ket {k 3 ,k 4 ) = Flip(A:i, ^ 2 ). 
Then 

\x{ki) - x{k2)\ = \x{k3) - x{k4)\- 

If, in addition, Z\,Z 2 ^ K and x(fci) x{k 2 ), then the following holds: 

\x{ki) - x{k2)\ yf \x{ki) - xiks)]- 
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Algorithm 2 FlipLinkage 

Input: ki,k 2 S K such that pp(fci) ^ pp(fc 2 ) and pp(fci), pp(fc 2 ) G C \ M. 
Output: {G,p) a flip linkage. 

1 ; Define V = {1,2,3,4} and E = {(1, 2), (2,3), (3,4), (1,4)} and set G = (V,E). 
2 : Define {k^,k 4 ) = Flip(fci,/c 2 ). 

3: Define Ui = xi^i) for alH G (1,..., 4}. 

4: Define 

p(l,3) = p(2,1) = f5i(u2), 9(u2)) 

p(3,4) = (3 ?(m3),9(u3)) p(4, 2) = (3?(u4),9(u4)) 

5: Return {G,p). 


Proof. Let ki = zi+pwi and ^2 = Z 2 + ?? W 2 • Then the condition (t—ki)(t — k 2 ) = 
{t — k^){t — kjf) implies that k^ = Z 2 + pws and A :4 = zi + 77^4 with 

- (-^2 -^)wi + {zi- Z^) W2 _ (Zl - Z 2 ) Wi + (Zi - Zj) W2 

zi-z^ zi-z^ 

A direct calculation shows that 


|X(^3) - X(fc4)| = 


W 3 

Z2 


Wi 

Zi- Zi 


Wi 

Zi- Zi 


W2 

Z2 


= \xiki) - Xik2)\ 


as claimed. The second claim is also obtained via a direct computation. 


□ 


Lemma 6.5 implies that the linkages produced by the algorithm FlipLinkage 
have the shape of a (possibly degenerated) antiparallelogram, as in the right part 
of Figure We introduce the notion of flip mobility to exclude some degenerated 
cases, namely the ones where the antiparallelogram is in fact a square. 


Notation. Let ki = zi + ijwi and k 2 = Z 2 + p 'W 2 be two elements of K. We say 
that the condition FM(fci, k 2 ) (which stands for flip mobility) holds if and only if 

Zi,Z2^K, ZX+Z2, Zxf^Z2, X(fcl ) 7^ xCfe ) • 

Lemma 6 . 6 . Lot fci,/c 2 G K such that FM(/ci, ^ 2 ) holds. Then the linkage L 
obtained by FlipLinkage(fci, fc 2 ) has mobility one and the configuration curve has 
two components Gi and G 2 . Moreover, the natural maps Gi —> VRP(i,j) and 
C 2 —> VRP(7, j) are isomorphisms for every two neighboring links i and j. 


Proof. Since the linkages obtained via FlipLinkage are constituted by four bars, 
this is a well-known result; see for example |31 Table on p. 426]. One can also prove 
the statement using the notions we introduced in Section]^ as follows. We set L = 
FlipLinkage(fci, A; 2 ), and denote (^ 3 ,^ 4 ) = Flip(A;i, ^ 2 ). Recalling Definition 3.6 
our goal is to prove that Conf(L) is a one-dimensional variety. One notices that 

Conf(L) = {(Ai :/ii),..., (A 4 : 774 ) : /i 2 (Ai,/ii, A 2 , 772 ) =/ 34 (A 3 , 7 x 3 , A 4 , 714 )}, 


where 

fi2: Pj^ X Pj^ > P(1K) 

(Ai : pi), (Ai : pi) (Ai — pi ki){X2 — P2k2) 

and analogously for f^^. The images of /12 and /34 are two smooth quadrics whose 
intersection, equal to Conf(L), is a curve with two components. Using this it is 
possible to prove the claim. □ 
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Remark 6.7. Results in Example |4.4| are particular instances of Lemma | 6 . 6 | One 
can in fact check that in that case L = FlipLinkage(fci, ^ 2)1 where ki = i — rj ■ Z 
and k 2 = —— ry (2 + ^). 


7 . Construction of linkages 

In this section, we show how to construct a linkage with mobility one that traces 
an algebraic curve described by a rational parametrization. Based on the results 
of Section we first describe an algorithm which takes a bounded motion polyno¬ 
mial P and constructs a linkage weakly realizing P (see Definition 4.141. 


Algorithm 3 ConstructWeakLinkage 

Input: P G K[t] a bounded motion polynomial. 

Output: {G,p) a linkage weakly realizing P. 

1 : Compute S G M[t] of maximal degree, monic, such that S divides P. 
2: Compute (R, {ki ,..., fc„)) = FactorMotionPolynomial(P/S'). 

3: Set Ui = x{ki) = ki = Zi + rjiVi, for every i S {1,..., n}. 

4: Set V = {1,..., n -F 1} and E = {{1, 2), (2,3),..., (n, n -F 1)}. 

5: Set G={V,E). 

6 : Set p{i, i -F 1) = (jR{ui), ^{ui)) G for every i S {1,..., n}. 

7: Return {G,p). 


Proposition 7.1. Algorithm ConstructWeakLinkage is correct. 

Proof. The linkage L constructed by the algorithm ConstructWeakLinkage is an 
open chain, so its configuration space is isomorphic to where n — 1 is the 

number of joints of L. Let (p: K —P^ be the motion corresponding to RP. We 
consider links 1 and n, and we prove that image((/)) C RP(l,n), namely that L 
weakly realizes RP according to Definition |4.14[ In our situation we have that 

RP(1, n) = {tTi ,2 O 0 - 2,3 O • • • O e Pc}. 

Now we fix an arbitrary t gM. and we take o-iy+i to be the isometry given hy t— ki 
for alH G {1,..., n — 1}. Then (j){t) G RP(1, n), so the claim is proven. □ 

Suppose we are given a rational parametrization of a bounded curve. As we saw 
in Proposition |4.7| starting from it we can construct a bounded polynomial P such 
that the curve is the orbit of the origin in under the motion described by P. Then 
the algorithm ConstructWeakLinkage returns a linkage L that weakly realizes P. 
However, this is rather unsatisfactory, since the configuration space of L is, in 
general, a variety of high dimension. Our goal is to obtain a linkage with mobility 
one, providing a strong realization of P. The main idea is to “rigidify” L, which by 
construction is an open chain, by introducing additional links and joints, forming 
antiparallelograms. At the level of graphs, this corresponds to extending the linear 
graph of L to a 2 X n ladder graph (see Figure]^. To do so, we pick a linear motion 
polynomial t — l, which will constitute the first “step” of the ladder, and then apply 
the flip procedure iteratively to complete the ladder. In the following, we investigate 
how I G K has to be chosen so that the resulting linkage has mobility one. 
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Definition 7.2. Let k = (fci,..., km) be a tuple in K™, and let I G K. Define 

li := I, (ki,li+i) := F\ip{k, ki) for alH G {1,..., m}. 

Then we say that the condition IFM(Z,k) (which stands for iterated flip mobility) 
holds if and only if FM{li, kf) holds for all i G {1,, to}. 

In order to prove that the algorithm we propose works, we have to show that 
it is always possible to find an element I G K such that the iterated flip mobility 
condition holds. For this we use a property of flip linkages. 

Definition 7.3. Given an element k = z + rjw G K with z yf 0, we define inv(A:) 
to be the element z — rj w G K. Recalling Remark |4.2| we have that k represents an 
isometry tr, and inv(fc) its inverse Moreover, if z ^ K, then t — k represents 

an isometry for every f G M, and t — inv(fc) the inverse isometry. 

Lemma 7.4. Let ki,k 2 G K be such that pp(fci) pp(fc 2 ) md pp(fc 2 ) ^ M. Let 
(k^jki) = Flip(fci,fc 2 ). Then we have that 

(inv(fc3),fci) = Flip(fc4,inv(fc2)). 

Proof. Notice that the fact that pp(fc 2 ) ^ M implies that the element {t — kf) 
represents an isometry for all t G K. From Definition |7.3| one has that the inverse 
isometry is represented by {t — inv(fc 2 )). The two situations are depicted below: 

t—fc2 t—inv(fc2) 


t—k4 t—ki 


t—k^ i—inv(fc3) 

[kz^kf) = Flip(fci,fc 2 ) (inv(fc3),/ci) = Flip(/c4, inv(fc 2 )) 

A direct computation shows the desired result. □ 


A 

t — k^ 


t — k-\ 


Lemma 7.5. Let k = (fci, ..., k„i) be a tuple in K'" such that pp(fci) ^ K for all 
i G {!,..., to}. Then there exists I G K with pp(Z) ^ M such that IFM(^, k) holds. 


Proof. Following the notation of Definition |7.2| we have that pp(^i) = pp(0 for all 
i G {!,..., to}. Hence, in order to ensure that IFM(/, k) holds, we need to choose I 
such that pp(/) y^ pp{ki) and pp(/) yf pp(fci) for all i G {!,... ,to}. From now on, 
we fix such a value for pp(Z). The other situations we should keep away from are 
the ones where xih) = xi^i)- Since pp(Zi) is now fixed, the condition x{h) = x(fci) 
becomes a linear equation for sp{li), which alw ays a dmits a unique solution uf G C. 
Let us define f — pp(0 + rjWi. Using Lemma 7.4 repeatedly (namely i — 1 times) 
starting from the pair (Z^, inv(fci_i)), we obtain an element G K such that 


x{h) 7 ^ xih) I h,. 

Summing up, for every fixed value of pp{l) different from all {pp(fci)}^j^ and all 
{pp(A;i)}™^, there are finitely many values for sp(Z) that should be avoided for 
IFM(Z,k) to hold. Thus we can always find I G K with the claimed properties. □ 


We are ready to present the algorithm ConstructStrongLinkage. Figure 
shows the labeled graph produced by the algorithm. 
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Algorithm 4 ConstructStrongLinkage 
Input: P G K[t] a bounded motion polynomial. 

Output: {G,p) a linkage strongly realizing the motion induced by P. 


Compute S G 


of maximal degree, monic, such that S divides P. 


Compute (i?, k = (fci,..., k^)) = FactorMotionPolynomial(P/S'). 
Choose I G K such that IFM(^, k) holds. 

Set Zi = Z and (fci, = Flip(Zi, Zci) for all z S {1,..., u}. 

Set 


Ui = Xih) 
Ui = Xih) 

vj = xih) 


for every i G {1,... ,n} 
for every j G {1, ..., n + 1} 


6: Set y = {!,...,2n +2}. 

7: Set B = |(z, z+1), (n+l+i, n+2+i), (j, n+l+j) for alH S {1,..., n} and j G 


8: Set G={V,E). 

9: Set 

p{ip + 1 ) 
p{n + l+z,n + 2 + z) 
p{j,n + l+j) 

10: Return (G,p). 




for every i G {1,... ,n} 
for every j G {1, ..., n + 1} 



Figure 7. The labeled graph of the linkage returned by the algo¬ 
rithm ConstructStrongLinkage. 


Theorem 7.6. Let P be a bounded polynomial in IK[t] and let 4>: M —> Pj| be the 
rational motion induced by P. Let us denote by L the linkage obtained by applying 
ConstructStrongLinkage to P. Then L strongly realizes the motion (j). 


Proof. Let K be the configuration space of L = {G,p). As in the statement of 
Lemma |6.6| we are going to prove that for any two neighboring edges i and j of G, 
every irreducible component of K is isomorphic to RP(z, j). This will prove that K 
has dimension one, which means that L has mobility one. 

Recall from Definition 3.4 that K is the collection of configurations E = 
satisfying the equations imposed by the directed cycles of G. First of all, notice 
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that one can consider only those equations imposed by directed cycles of the form: 

h • -s- 9 h + 1 


{n + l) + h 9^ -• (n + l) + (/i + l) 


In fact, all other equations belong to the ideal generated by these ones. We intro¬ 
duce the following notation: for h G n}, denote by sq(/i) the set 

We define Trsq(h) to be the projection 


T^sq{h) ■■ K —!■ VRP(r, s). 

(r‘,s)Gsq(/i) 


Then the image o f 7rsq (^) is contained in the configuration space Ch of a flip linkage, 
which by Lemma 6.6 has two components, so Ch = C^- 

For every sequence 6 S {1, 2}”, we set 


Kb 



By construction we have K = lj&6{i 2 }^ We are going to prove that for every 
sequence b and for any two neighboring links i and j in G, we have Kb = VRP(*, j). 

For the sake of simplicity, we start by proving this result when the graph G 
consists only of two squares (i.e. n = 2): 



In this case we only have two flip linkages with configuration spaces Gi and G 2 . 
Let us fix 6 S {1, 2}^, and consider the diagram: 



VRP(2,5) 


The maps 712,5 ar® isomorphisms, because of Lemma 6.6 We consider the set 
|(a:,y) : a: € Gi^^\y S G^^^^ and 712 , 5 ( 0 ;) = 712 , 5 (j/)|, 


which is nothing but the pullback of Diagram One can check that this set 

equals Kb since, as we observed before, the elements of K are only subject to the 
equations coming from the cycles in Gi and G 2 . Since isomorphisms are stable 
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under pullbacks, we obtain the following commutative diagram of isomorphisms: 

Kb 



VRP(2,5) 


The composition 712,5 o Tsq(i) equals the projection 712 , 5 : Kb —> VRP(2, 5), and 
the same holds for every projection T^ij: Kb —>■ VRP(i,j), where i and j are 
neighboring links. So each of these maps is an isomorphism, proving our claim. 

If now G is constituted by more than two squares, then Kb is obtained via an 
iteration of several pullbacks, and eventually we get a diagram of isomorphisms. 
Thus the claim holds also in this case. Since K is covered by finitely many varieties 
of dimension one, it has dimension one. Moreover, each of the varieties Kb is an 
irreducible component of K. 

We are left to prove that L realizes the motion (j) induced by P. Similarly as we 
did in the proof of Proposition [Z3 one can show that image{(/>) C RP(l,n + 1). 
Since K is one-dimensional, also RP(l,n-|- 1) is so; since image((/)) is one-dimen¬ 
sional, image(0) is a component of RP(l,n -|- 1). This concludes the proof. □ 


Remark 7.7. As we mentioned in Section the linkages constructed via our 
algorithms present the same issues as the ones produced by Kempe’s procedure, 
namely the devices do not realize only the motion they are designed for, but also 
other ones. This is due to the fact that the elementary linkages produced by 
the flip procedure, commonly known as antiparallelograms, admit configuration 
spaces with more than one component (see Lemma 6 . 61 . On the other hand, since 
antiparallelograms are among the linkages used also by Kempe’s procedure, the 
techniques developed in |T] and O Section 3.2.2] can be applied also in our case to 
prevent this unwanted behavior. 


Remark 7.8. Let (j)he & rational motion given hy Z + rjW G IK[t]. Our algorithm 
allows one to construct a linkage that generates (j), using a motor that rotates at 
constant speed. For this purpose one needs two neighboring links whose relative 
position is a revolution with constant speed; such a motion is represented by a 
motion polynomial of the form t ± i + rjw. If z or —j is already a root of Z, then 
any linkage returned by ConstructStrongLinkage has the desired property. If not, 
then the linear polynomial I in step 3 of this algorithm has to be chosen accordingly. 


We now investigate how many links and joints the linkages created by our al¬ 
gorithm have. We consider the tasks of strongly realizing a given rational motion, 
and of drawing a given parametrized curve. 


Proposition 7.9. Let (j) be a rational motion given by a motion polynomial P of 
degree d. Then there exists a linkage with at most Ad + 2 links and 6d -I-1 strongly 
realizing (j). 

Proof. Apply ConstructStrongLinkage to P, and notice that the real polyno¬ 
mial R produced by FactorMotionPolynomial has degree at most d. □ 
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Proposition 7.10. Let (p: M — > be the parametrization of a real bounded 

planar curve. Without loss of generality, suppose that 

9if)\ 

\h{t)’ h{t)) 

with f,g,h G and h monic, and d := degh > inax{deg/,deg(/}. There exists 
a linkage with at most 3d + 2 links and | d + 1 joints drawing the curve given by (p. 


Proof. As discussed before Definition 5.13 the motion polynomial P = h+p {f+ig) 
is bounded, and by Proposition 4.7 the orbit of the origin in under the motion 
determined by P is the image of ip. We could apply ConstructStrongLinkage 
to P, but we obtain a smaller linkage by applying this algorithm to P' — CP, 
where C G C[t] is any polynomial such that CC = h and gcd(C, C) = 1; such a C 
exists because h has no real roots. In fact, as discussed in Remark |4.9[ the orbit 
of the origin under the motion given by P' is also the image of ip. Moreover, by 
Theorem |5.15| a motion polynomial Z -\- pW is factorizable if the corresponding 
real polynomial R equals 1. This is the case if and only if gcd{Z,Z) \ WW; it 
is easy to check that this condition holds for P'. Finally, the linkage obtained by 
applying ConstructStrongLinkage to P' is constituted by 3d + 2 links and | d + 1 


joints, because the degree of P' is | d 


note that d must be even. 


□ 


Example 7.11. We illustrate the algorithm ConstructStrongLinkage with the 
example from Section In contrast to Example |5.16| here we multiply the poly¬ 
nomial P(t) = (t^ + 1) + p (it — 2) by C = t — i, as prescribed by Proposition 


7.10 


Calling FactorMotionPolynomial with CP as input, we get i? = 1 and obtain: 

C{t) ■ P{t) = (t — i — ^pi) ■ (t — i + ^pi) ■ (^t + i + pi). 

Next we have to fix an / G IK such that the IFM condition holds; it turns out that 


I = — ^i — T^p i = li is a good choice (which yields a linkage that is well suited for 


9 „ 18 „ 

/ t — 

visualization, see Figure]^. Applying the flip procedure iteratively to this input 
data, one obtains the following: 


ki = i — ^pi, 


k2 = ' 




l2 = -p 


“cl — 5“ 40'/ “) “4 — 5 . 

Now we can use these quantities to construct a linkage — whose graph is shown 
in Figure — that draws an ellipse. Note that all linear polynomials appearing 
here are purely imaginary in their secondary parts. This implies that their fixed 
points are of the form (cc, 0), so that all joints are located on the horizontal axis 
when the linkage is in its initial position (t = 00 ); the same happens for t = 0. 


h = -h- 
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k3 = -i- %pi 


Ia = -\ 


fi^*- 



Figure 8. The link graph for the linkage drawing an ellipse. 
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8. Self-collisions 


A natural question that arises in the construction of a linkage is whether it can 
be physically realized such that no collisions between components of the linkage 
occur. More precisely, one asks whether there exists an assignment of 1,2, ...,n 
to the links such that for any i < k < j for which i and j are neighboring links, 
the joint connecting i and j never overlaps with the link k. For general linkages it 
is a difficult problem to detect such collisions; see for example [51 Section 9.3 and 
Theorem 9.5.5]. However, it turns out that the same problem has a straightforward 
solution for linkages constructed by our algorithm, and the reason is that all joints in 
our case follow rational curves whose explicit parametrization is a direct byproduct 
of the algorithm ConstructStrongLinkage. Assuming that a link is realized as 
a collection of line segments connecting the (two or three) joints attached to it, 
a collision is described as follows: for some i < k < j as above there exists a 
t e K U {oo} such that the position (^xi{t),yi{t)) of the joint {i,j) lies on one of 
the line segments of link k; denote its endpoints by [x 2 {t), y 2 {t)) and {x 3 {t), ysit)). 
Then a collision between joint {i,j} and link k happens if and only if the system 


Xi{t) = S ■ X2{t) -f (1 - s) • X3{t) 
yi{t) = s ■ y 2 (t) -k (1 - s) • y^it) 


has a solution for 0 < s < 1 and t G M.U {oo}. In our construction, the coordi¬ 
nates Xi,yi are rational functions in t, so that the system ( 8.11 is equivalent to a 


bivariate polynomial system that can be easily solved. Note that we have to solve 


(8.1 1 for each admissible triple i < k < j. This way we can decide whether collisions 


occur or not, and if so, we get a precise description when and where they happen. 


Example 8.1. Consider again the linkage constructed in Example 7.11] whose link 
graph is shown in Figure With the ordering (5,1, 6 ,2, 7, 8 ,4, 3) of the links, we 
get two collisions, which both happen at t = oo. Hence in principle we can trace, 
without disassembling the linkage, the full ellipse except a single point; see also the 
animations provided on our webpage HZ]. 


At the same time, we want to mention a remarkable method for preventing 
collisions for planar linkages with a ladder-shaped link diagram as shown in Figure]^ 
This method is based on designing the shape of some links and joints: we introduce 
three types of links (F-link, U-link, Z-link), and two types of joints (T-joint, Z-joint). 

As before, all links are arranged in different layers, and w.l.o.g. we assume that 
these layers correspond to integer numbers. An F-link is located in a single layer, 
in other words, we associate one integer for each F-link. On the other hand, U-links 
and Z-links stretch across two, not necessarily neighboring, layers. Therefore, to 
such links we associate a pair of integers (a, b) with a < b, which means that one 
part of the link is located in layer a and the other part in layer 6 ; the two parts are 
rigidly connected by a vertical rod. For two U-links with layers (oi, 6 i) ^ ( 021 ^ 2)1 
we prohibit the situations oi < 02 < 61 < 62 and 02 < oi < 62 < &i, which could 
yield a collision between these two U-links. For a Z-link on layers (a, b) we impose 
the condition b — a = 2, and the link located on layer a -I- 1 (this will always be 
some part of a U-link) is connected with this Z-link by a revolute joint (which we 
call Z-joint) around its vertical rod. In contrast, a T-joint joins two links located 
on neighboring layers. These types of links and joints are illustrated in Figure]^ 
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Using the above design, one can check that collisions can only happen between 
links, not between links and joints. We argue now that even such collisions can be 
avoided. First, an F-link cannot collide with another F-link or a Z-link. Second, 
a Z-link cannot collide with another Z-link. Third, if we move the vertical rods of 
all U-links sufficiently far away, they do not collide with the F-links and Z-links. 
Fourth, the above conditions on the layers of U-links imply that two U-links can 
collide only if they occur in a nested way (e.g., Oi < 02 < &2 < ^i)- Again, this can 
be avoided by moving the vertical rod of the outer U-link far enough away. Thus we 
can always make a collision-free design, by manipulating the shape of the U-links. 



Figure 9. A linkage without self-collisions composed of three F- 
links (gray), two U-links (white), and one Z-link (black) 


It remains to argue that for a general ladder-shaped link diagram, as returned 
by our algorithm ConstructStrongLinkage and as displayed in Figure we can 
find an assignment of link types and layers such that all the above conditions are 
fulfilled. One such assignment is depicted in Figure it demonstrates that we 
can realize a motion — given by a factored motion polynomial of degree n — by a 
collision-free linkage using 4n -|- 1 layers. 


(0, 4n) 


(3,4n-l) (6,4n-2) 


(3n—3, 3n-|-l) 


(3n) 



( 1 ) 


(2,4) 


(5,7) 


(3n-4,3n-2) (3n-l) 


Figure 10. Assignment of joint types (T, Z), link types (F, U, Z) 
and layers (above and below the corresponding vertices) for realiz¬ 
ing a linkage with ladder-shaped link graph without self-collisions. 
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Figure 11. A rational curve approximating the “J” in John Han¬ 
cock’s signature and a linkage drawing it. The corresponding mo¬ 
tion is a translation, which we visualize by using a quill pen whose 
shape is a line segment in direction (5,6). The spatial arrangement 
of the links is indicated by hue: darker links lie below brighter ones. 

9. A SIGNING LINKAGE 

We conclude with an example in connection to a popular formulation of Kempe’s 
Theorem, stating that “There is a linkage that signs your name”. King [T?l Corol¬ 
lary 1.3] attributes this formulation to William Thurston. However, as remarked 
by O’Rourke [53], it is very implausible that a concrete “signing linkage” has ever 
been constructed due to the complexity, in terms of links and joints, of the link¬ 
ages produced following Kempe’s procedure. As an example to support his claim, 
O’Rourke points out that already constructing a linkage drawing the “J” of John 
Hancock’s famous signature (see Figure [T^ would be very difficult. 

We approximate the “J” by the rational curve given by the parametrization 
{flh,glh), where 

f{t) = -321880t'^ - 436132t^ - 237449^^ - 6U88t^ - 8666t - 451, 
g{t) = -336018t® - 472949t^ - 270569t^ - 78158t^ - 11325t - 651, 
h{t) = 170 (7225t® -k 13770t^ -k 11187t'‘ -k 4908f^ -k 1219t^ -k 162t -k 9) . 

Our implementation m of the algorithm ConstructStrongLinkage produces a 
linkage with 26 links and 37 joints realizing a translational motion along this curve, 
see Figure mi notice that these numbers are in accordance with Proposition |7.9| 
When this linkage draws the depicted “J”, one encounters seven collisions. 
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Figure 12. John Hancock’s signature on the United States Dec¬ 
laration of Independence. 
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